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ABSTRACT 

We revisit the apparent correlation between soft X-ray band photon index and 
spin-down rate i> previously report ed for Anomalous X-ray P ulsars (AXPs) and 



Soft Gamma Repeaters (SGRs) by iMarsden White! (120011 ). Our analysis, im- 



proved thanks to new source discoveries, better spectral parameter measurements 
in previously known sources, and the requirement of source quiescence for pa- 
rameter inclusion, shows evidence for the previously noted trend, although with 
greater scatter. This trend supports the twisted magnetosphere model of magne- 
tars although the scatter suggests that factors other than i> are also important. 
We also note possible correlations involving the spectra of AXPs and SGRs in the 
hard X-ray band. Specifically, the hard-band photon index shows a possible cor- 
relation with inferred v and B, as does the degree of spectral turnover. If the for- 
mer trend is correct, then the hard-band photon index for AXP IE 1048.1—5937 
should be ~0-l. This may be testable with long integrations by INTEGRAL, or 
by the upcoming focussing hard X-ray mission NuSTAR. 

Subject headings: stars: neutron — pulsars: general — X-rays: stars 



1. Introduction 

Anomalous X-ray Pulsars (AXPs) were long thought of as "soft" X-ray sources, because 
their X-ray spectra below 10 keV were observed to be falling off sharply. INTEGRAL 
and Rossi X-ray Timing Explorer (RXTE) observations of these objects, as well as of their 
close cousins, the Soft Gamma Repeaters (SGRs), revealed that their spectra turn over in 
the 10-20 keV range, such that in fact some AXPs are the hardest known sources in the 
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sky in the hard X-ray band (|Molkov et al.l l2004 iKuiper et al.ll2004 iMereghetti et al.l 12005 
den Hartog et al.ll2006t IGotz et al.ll2006l ). Here we define the "soft" X-ray band as 1-10 keV, 
and the "hard" X-ray band as 20-80 keV. 



Prior to the discovery of the hard X-ray emission, iMarsden fc White! (120011 ) (hereafter 
MW01) reported an apparent correlation between soft-band spectral hardness and spin-down 
rate z>, where v is the spin frequency. This reported correlation has b een used as the primary 



motivation for the "twisted magnetosphere" model for magnetar s ([Thompson et al.l 12002 



Thompson fc Beloborodovl l2005t iBeloborodov fc Thompson! 120071 ) . In that model, surface 
activity driven by the decay of the large internal magnetic field of the star results in a twisted 
dipolar morphology of the magnetosphere, with powerful currents emanating from the sur- 
face and subsequently returning to heat it. In this picture, the MW01 hardness/z> correlation 
is interpreted physically as a correlation between hardness and spin-down-inferred surface 
dipolar magnetic field B, where the latter can be estimated via B = 3.2 x 10 19 (PP) 1 ' 2 G, 
with P = 1/u the pulse period. The simple z>/B connection exists because, for AXPs and 
SGRs, there is a far wider range in observed v than observed v. In the twisted magneto- 
sphere model, this correlation is a result of stronger currents, due to higher B hence higher 
activity, resonantly upscattering seed thermal photons from the hot surface to form the non- 
thermal component. In addition to the apparent soft-band hardness/S correlation reported 
by MW01, support for this picture comes from observed correlations between soft-band spec- 



tral index and flux in individual varia ble sources (e.g. iTiengo et al.l 120051 : iRea et al.l 12005 



Campana et al.ll2007l ; iTam et al.l 120081 ). since all else being equal, higher fluxes are expected 
from larger twists, because of higher return currents. 

In this Letter, we revisit the correlation reported by MW01 in light of three main de- 
velopments in the field since their work. First, several new sources have been discovered 
since 2001 and can be used to confirm the original claimed correlation, which was based on 
only seven objects. Second, previously known sources have been observed in longer inte- 
grations with higher sensitivity instruments, allowing better determinations of their spectral 
parameters. Third, since 2001, it has been recognized that AXPs can exhibi t great vari- 
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measured magnetar spectral parameters without concern for the objects' outburst status at 
the time of measurement, thereby possibly using data from objects in very different spectral 
states. 

In addition, with the study of the hard-band emission from AXPs and SGRs well un- 
derway, there now exist enough data for correlations similar to that noted by MW01 to be 
considered. Indeed IGotz et al.l (120061 ) noted qualitatively that in the hard band, AXPs are 
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harder than SGRs, the reverse of what is seen in the soft band. Here we consider this point 
in more detail, noting possible correlations between the hardness of the hard-band spectra 
of AXPs and SGRs and B (and u), and between magnitude of spectral change from the soft 
to hard bands, and B (and v). 



2. Compiling AXP and SGR Spectra 



We characterize AXP and SGR soft and hard X-ray spectra phenomenologically. Al- 
though signific ant effort has been expended to do physica l ly motivated broadband spectral 



modelling fe.g. iLvutikov &: Gavriilll2006t iRea et al 



2007b: Ozel & Guver 2007; Nobili et al. 



2008 ; Zane et al. 20091 ). this work is still progressing, with no final word reached. Moreover, 



our motivation here is to characterize all sources in a systematic and uniform way. Hence we 
choose to parametrize soft-band spectra using the simple blackbody plus power-law model, 
with T s being the soft-band photon index (defined as dN/dE oc E~ Ts ). Similarly, we param- 
eterize the hard-band spectrum by a power law of photon index (defined as for T s ). We 
recognize that this parameterization is flawed, both in the soft band where the true spectrum 
is likely to be better represented by a Comptonized thermal model (such as those currently 
under development) as required in the twisted magnetosphere model, but also more dramat- 
ically in the mid-energy range where there is clearly spectral curvature. Nevertheless the 
parameterization is convenient, used ubiquitously in the literature, and does a reasonable 
job of characterizing the soft and hard-b and spectra currently available, given statistical 
uncertai nties (see e.g. iKuiper et al" 2006 ) . We further acknowledg e that the spectra, both 



soft (e.g. iPatel et al.ll200ll ) and hard (e.g. Iden Hartog et al.ll2008al ). are dependent on pulse 



phase, however here we consider only phase-averaged spectral parameters. 

In Table 1, we summarize the results reported in the literature for all known non- 
transient magnetars for which there are such measurements. Since MW01, several new mag- 
netars have been discovered, four of which have well measured spin and spectral parameters 
(CXOU J010043. 1-721134, SGR 0526-66, IE 1547.0-5408 and CXO J164710.2-455216). 
In addition, in the soft X-ray band, SGRs and AXPs are now established to be highly vari- 
able, such that T s can vary greatly. In compiling Table 1, we took great care to include the 
T s (always for a blackbody plus power-law model) for objects when their reported fluxes 
were at or near the lowest yet recorded. This was done to ensure the objects were in their 
quiescent state, with no temporary enhancement that could be interpreted as a magneto- 
spheric twist with corresponding spectral hardening. By contrast, MW01 averaged together 
different published values of T s , presumably under the assumption that the values do not 
change, which we now know to be incorrect. In making our selections from the literature, 
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where multiple quiescent spectral observations were available, we verified that the reported 
T s values were roughly consistent, then selected that with the smallest reported uncertainty. 
We excluded the so-called transi ent AXP XTE J1810 — 197 which shows no clear non-thermal 
component when in quiescence ( IGottheff et al.ll2004f ). 



3. Soft-Band Spectral Indexes 

Using the data provided in Table 1, we plot T s versus B and i> in Figure [TJ Values 
of v are long-term averages as recorded in the SGR/AXP online catalogs! (see references 
for z> values therein). We note the previously seen trend for decreasing T s with increasing 
B or v. To quantify this trend, we determined Pearson's linear correlation coefficient to 
be r = —0.56 for a sample size of iV = 11, corresponding to a null-hypothesis (two-tailed) 
probability of p = 0.076. We note that this p value, as for all others in this paper, is actually 
an upper limit given that uncertainties are not accounted for in calculating r. This T s /z> 
trend is not as striking as it was in MW01 (compare with their Fig. 2); there is more 
scatteJE If we ignore AXP IE 1547.0—5408, whose T s has large uncertainty, the trend is 
improved and is significant with > 99% confidence (r = —0.78 for N = 10, p = 0.0079), 
though still with considerable scatter. Interestingly, the plot of T s versus B (Fig. (TJ left) 
shows a better correlation (r = —0.82 for N — 11, p — 0.0022), which may support the 
twisted magnetosphere model in that it is the magnitude of B that determines the twist 
size, not just the spin-down rate. Nevertheless, there is scatter; this could imply that the 
value of B, as inferred from spin-down parameters, is only a rough estimate of the true field 
strength, a possibility already generally recognized given the greatly simplified assumptions 
in computing B, and also observationally suggested from the absence of anomalous X-ray 
emission from very high-5 rotation-powered pulsars (e.g. Pivavoroff, Kaspi & Camilo 2000; 
McLaughlin et al. 2003; Kaspi & McLaughlin 2005). 

We also searched for a correlation between the simultaneously measured blackbody 
temperature kT and B (see Fig. 2), but, like MW01, found none. Specifically, for N — 11, 
r = 0.29, implying p = 0.38, clearly non-excluding of the null hypothesis. This argues that 
the underlying physical parameter that is most closely linked with a magnetar's magnetic 
field is manifested in the power-law component, rather than in the thermal component (see 



also iKulkarni et all 120031 ) . 



1 www. physics. mcgill.ca/^pulsar/magnetar/main. html 

2 We estimate for MWOl's reported T s vs. v correlation (their Fig. 2), r = —0.83 for N — 7, implying 
p = 0.021. 
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4. Hard-Band Spectral Indexes 

Figure [3] shows values of i\ plotted versus B (left) and v (right) for the six sources for 
which Th has presently been measured. In some cases, Th is for the pulsed emission, in other 
cases for the total emission, and in three cases for both (see Table 1). There is a possible 
trend in the data. AXPs ar e harder in the 20-80 keV band than are SGRs as qualitatively 



noted by iGotz et al.l ( 120061 ) . But as is clear here, even within the AXP class itself, low B 
or v sources are harder than high B or v sources. To be more quantitative, for the sources 
for which the pulsed emission is measured, we find r = 0.88 for N = 4, implying p — 0.12 
for a r\ correlation with logB (Fig. 3, left), and r = 0.86 for N = 4 implying p = 0.14 for 
the T/j versus logz> data (Fig. 3, right). These values are more significant for the total flux 
(r = 0.79, N = 6,p = 0.062 and r = 0.84, N = 6,p = 0.034 for log.B and logz>, respectively) 
when the (pulsed) value for IE 2259+586 is included, under the assumption that its total 
flux Th does not differ much from its pulsed Th, as for the three sources for which both 
are measured. Although not quite significant at the > 99% confidence level, these possible 
correlations seem noteworthy. 

We note that the value of Th for IE 2259+586 should be regarded with caution, the 
source having been detected only up to 30 keV thus far. This will hopefully improve with 
time. We also of course caution that the trend is defined by only six sources, one fewer 
than was used by MW01. If this correlation is correct, then for AXP IE 1048.1—5937, the 
brightest AXP for which Th has not yet been measured, Th should be ~0-l. This may be 
testable with long integrations by INTEGRAL, or by the future focussing hard X-ray mission 
NuSTAR (expected launch 2011) which should have point-source sensitivity of ~100 times 
that of INTEGRAL. 

Note that in all but two cases (4U 0142+61 and SGR 1806—20), the sources were in qui- 
escence at all epochs when the hard X-ray observations were done. For 4U 142+61, the mea- 



sured hard X-ray spectral parameters are an average over ~4 yr of data (Iden Hartog et al. 



2008bl ) which included a few-month active period in 2006 (jGavriil et al.l 120091 ) . This active 



period seems unlikely to have biase d the hard X-ray sp ectrum because of the absence of soft 



X-ray variability during that epoch (jGavriil et al.ll2009i ). because of the relatively small dura- 



tion of the hard X-ray exposure during the active period relative to during quiescence, as well 
as because of the absence of any evidence for hard X-ray flux or spectral va riability over the 4 



years o f INTEGRAL observations (with flux and stable to within 17%; Iden Hartog et al. 



2008b|). For SGR 1806—20, the hard X-ray observations were obtained roughly 2 yr after 
the source's giant flare of 2004 December. During the hard X-ray observations, the source 
was nearly back to quiesence, though still bursting occasionally. Whether this mild activity 
affected its hard X-ray spectrum is unknown given that no variability has yet been reported 
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in the hard X-ray spectrum of any magnetar. This will likely be difficult to determine with 
current hard X-ray telescopes given their sensitivity, barring another major outburst. 

Figure H] shows values of r s — 1\, our parameterization of the degree of spectral turnover, 
plotted versus B (left) and i> (right). A trend is suggested by both plots, such that the lowest 
B or v sources generally have the largest spectral turnover. Quantitatively, for the pulsed 
data, we find r = —0.99 for N = 4 implying p = 0.014 and r = —0.98 for N = 4 implying 
p = 0.021, for T s — T h versus B and z>, respectively. When considering total flux, these 
values change to r = —0.90 for N = 6 implying p = 0.015 and r = —0.92 for iV = 6 implying 
p = 0.0087, for T s — Th versus B and z>, respectively, again including the (pulsed) data for 
IE 2259+586, under the reasonable assumption that the difference between its total and 
pulsed flux spectral turnovers is small (as it is for the three sources for which both have 
been measured). Thus, these correlations are significant at the £ 99% level. If the versus 
B trend holds and the above prediction for Th for IE 1048.1—5937 is correct, its spectral 
turnover would agree well with this second trend. The same obvious caution regarding the 
paucity of sources as in Figure [3j holds, however. 

We also investigated whether there exists any correlation between hard X-ray luminosity 
and B, but found none. However, the luminosities of the sources in this band are generally 
poorly known, given relatively large Th and distance uncertainties. Thus some correlation 
may exist, but cannot presently be discerned. 



Discussion 



We have revisited the correlation between soft X-ray band power-law spectral index and 
i> and B reported previously by MW01, and find it continues to hold, though with significant 
scatter. This stands in contrast to the absence of any apparent relationship between black- 
body te mperature and B or i>. T his supports the foundatio n for the "twisted magnetosphere" 
model of T hompson et al.l (120021 ) , as fu rther developed by lThompson &: Beloborodovl (120051 ) 
and iBeloborodov &: Thompson! (120071 ). The significant scatter in the still small measured 
population suggests in the context of this model that the inferred B is only a rough estimate 
of the stellar dipolar magnetic field, and/or that other factors, such as geometry (both intrin- 
sic and viewing) and higher-order magnetic multipoles, play physically importan t roles. With 



very recent discoveries of se yeral new magnetars, for e xample SGR 0501+4516 (jEnoto et al. 



2009h and SGR 0418+5729 Ivan der Horst et al.lbood ). as well as new measu rements of spin 



param eters for previously known objects, as, for example, for SGR 1627—41 (jEsposito et al. 



20091 ). it should be possible to further test the veracity of the correlation in the near future. 
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The possible correlations we have noted for the hard-band emission from magnetars 
may be telling us something new about the magnetar emission process. Next we consider 
whether the published models have predicted the observed trends. 



Heyl fc Hernquistl (120051 ) suggested that magnetohydrodynamic fast modes produced 
by magnetars may break down to form electron/positron pairs in the magnetosphere. They 
suggested that non-thermal radiation associated with this fast-mode breakdown could ac- 
count for the high-energy emission from magnetars. However in this model, the spectrum 
of the radiation is fixed and is not expected to var y from source to source as has been ob- 
served. This difficulty, along with the now falsified (Iden Hartog et al.ll2008bl ) prediction that 
the fast-mode breakdown radiation spectrum should extend above 1 MeV, argues against 
this model being relevant to explain the observed properties of the hard X-ray emission in 
magnetars. 



Baring fc Harding! (120071 ) explore resonant inverse Compton upscattering of thermal 



emission from the surface, by relativistic electrons only a few stellar radii from the surface. 
This resonant scattering, at least in their preliminary analysis, is expec ted to yield a very 
flat sp ectrum, independent of B, unlike what is reported here. However, iBaring fc Harding 
(hoOTT l did not consider non-resonant scattering which they suggest may have a very different 
spectrum, although its relative contribution is as yet unclear. For the resonant scattering, 
although substantial spectral differences are expected for different scattering locales, as well 
as for more complicated field geometries than the assumed dipole, the flatness of the spectrum 
independent of B seems robust. The correlations we report thus suggest that this process is 
not the dominant mechanism for producing the hard X-ray emission in magnetars. 



hompson et al.l ( 120021 ) , iThompson fc Beloborodovl (120051 ) and lBeloborodov fc Thompson 



(120071 ) argue that powerful magnetic currents in the magnetar magnetosphere are induced 
by sporadic starquakes that twist the external magnetic field. A self-induction electric field 
lifts particles from the surface, accelerates them, and yields avalanches of pair creation. This 
forms a corona close to the neu tron star surface. IThompson &: Beloborodovl ( 120051 ) and 
Beloborodov &: Thompson! (120071 ) consider a model in which the hard X-rays are produced 
in the transition layer between the corona and thermal photosphere. In this model, the 
hard X-ray emission is a result of thermal bremsstrahlung with a temperature that scales 
approximately as kT oc B 2 ^ 5 , with kT ~ 100 keV for reasonable parameters. However, we 
verified, using spectral simulations, that this scaling predicts the opposite of the correlation 
seen in Figure |3j because higher kT, expected for higher B, gets fit with a lower when 
fit with a simple power law, for kT in the range 20-200 keV. Thus the observed correlation 
does not seem consistent with a thermal bremsstrahlung origin. 



On the other hand, the observed correlation between the spectral index change and B 
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could be interpreted at least qualitatively in the framework of models in which magneto- 
spheric currents play a role in the production of the high-energy emission as follows. In the 
magnetosphere, relativistic electrons can resonantly scatter surface thermal photons to form 



the no n-thermal portion of the soft-band spectrum as discussed first by iThompson et al. 



( 120021 ) . Alternatively, electrons that do not have a chance to sca tter will impact on the sur- 
face, h eating it, and possibly creating a corona as described by IThompson fc Beloborodovl 
(2005) and further explored by lBeloborodov fc Thompson! (120071 ) (but see lLyubarsky fc Eichler 
20071 ). For higher B so urces, the scattering opt ical depth increases, rendering the soft index 
harder as described by IThompson et al.l (120021 ). but leaving fewer electrons for heating the 
corona, and hence possibly resulting in the hard index being softer. 

For energetic rotation-powered pulsars (RPPs), the power-law spectrum seen below 
10 keV extends well beyond that energy, such that the spectral turnover r s — w for 
practically all sources for which it has been measured. In a handful of cases, there is marginal 
evidence for a subtle change of shape of the hard spectrum, gene rally parameterized as a 
log-parabolic steepening above 10 keV (e.g. ICampana et al.l 120081 ). Clearly if the trend in 
FigureHJheld for RPPs, one would expect huge turnovers, given their far lower B. This argues 
strongly that the mechanism for the production of hard X-rays in magnetars is different in 
RPPs, which was already suggested by the fact that the magnetar hard X-ray emission has 
luminosity far in excess of that available from rotational slow-down, in c ontrast to that from 
RPPs. T h for the usually rotation-powered PSR J1846-0258 (but see iGavriil et aJhoOsh. 
thehighest-S (4.9 x 10 13 G) RPP for which has been measured ( IKuiper &: Hermsen 
20091 ). is 1.2, much larger than would be predicted by the trend in Figure 3, left (T h versus 
B) but smaller than what would be predicted in Figure 3, right (Th versus v). Indeed the 
same is true of practically all RPPs for which Th has been measured, further demonstrating 
a different hard- X-ray production mechanism. 

We thank A. Beloborodov, D. Eichler, M. Livingstone, and C. Thompson for useful 
conversations. We thank Caltech for its hospitality and support via a Summer Undergrad- 
uate Research Fellowship (KB) and a Moore Scholarship (VMK). VMK receives additional 
support from NSERC, FQRNT, CIFAR, and holds a Canada Research Chair and the Lome 
Trottier Chair in Astrophysics and Cosmology. 
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Table 1: Spectral Parameters for Non-transient SGRs and AXPs 3 



Source Name 


B 


V 


BB kT 


r s 


Ref. 


K 


Ref. 


r* 

1 h 


Ref. 




(10 14 G) 


(IO" 12 s-2) 


(keV) 




for F 3 , kT 




for i* 


for r* 


IE 2259+586 


0.59 


-0.0099 


0.411(4) 


4.10(3) 


1 


-i.02g:ig 


12 






CXO J164710.2-455216 


~0.9 b 


-0.007 b 


0.49(10) c 


q k1.3c 


2 










4U 0142+61 


1.3 


-0.026 


0.386(5) c 


3.67(9) c 


3 


0.40(15) 


13 


0.93(6) 


13 


IE 1547.0-5408 


2.2 


-5.41 


n 4q0.03c 
U - 4C '0.04 


q 7 0.8c 
2.0 


•1 










CXOU J010043. 1-721134 


3.9 


-0.29 


0.35(6) c 


1.8(l) c 


5 










IE 1048.1-5937 


4.2 


-0.65 


0.56(1) 


3.14(11) 


6 










1RXS J170849.0-400910 


4.7 


-0.16 




9 7Q90.008c 
z - ,az 0.012 


7 


1.01(12) 


12 


1.44(45) 


12 


SGR 1900+14 


6.4 


-2.91 


0.47(2) c 


1.9(l) c 


8 






3.1(5) 


14 


IE 1841-045 


7.1 


-0.30 


0.44(2) 


2.0(3) 


9 


0.72(15) 


12 


1.32(11) 


12 


SGR 0526-66 


7.3 


-1.00 


0.48(5) 


3.12(8) 


10 










SGR 1806-20 


21 


-9.62 


0.476(39) d 


1.67(15) 


11 






2.0(2) c 


15 



Refs: [11 Woods et al. (2004): [2l Israel et aJ 


. (2007): [3] White et al. (19961): [41 Gelfand & Gaenslerl d2007l): 


[51 McGarrvetal. (20051): [61 Tarn et al. 




2008): [71 Rea et all d2007al): [81 Mereehetti et al. J2OO6I): [91 


Morii et al. (2003): [101 Kulkarni et al. 


20031): [111 


Woods et alJ ( 


2007); [12] KuineretalJ (20061): [13] 


den HartoE et al. (2008bl): [14] Gotz et al. 


(2006): [151 


EsDOsito et al. 


(20071) 



^Numbers in parentheses represent la uncertainties in the least significant digit(s) unless otherwise indicated 
(see c). 



h B and v for this source are not yet measured precisely; see lWoods et al.l ((2009). 
c These uncertainties represent 90% confidence intervals. 

d Thc BB temperature is from a different observation as the reported T s ; see lWoods et al.1 (|2007[ ). 



13 



3 



0526 J 



1 84 1 1 
1900 17 



-0.5 0.5 1 

Log(B 10 14 G) 



1.5 



3 



1048 J J 0526 



181 i 
0100 



-1 1 

Log(-i> 10- 12 s- 2 ) 



Fig. 1. — Left: T s versus B. See text and Table [5] for data used and references. Right: 
Same but versus i>. Both: Error bars represent la uncertainties. 
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Fig. 2. — Blackbody temperature versus B, for all magnetars for which T s is measured. See 
text and Table 1 for data used and references. Error bars represent la uncertainties. 
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Fig. 3. — Left: I\ versus B, for all magnetars for which I\ is measured. See text and 
Table 1 for data used and references. Circles represent pulsed flux, triangles represent total 
flux. Right: Same but versus z>. Both: Error bars represent la uncertainties. 
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Fig. 4. — Left: Spectral turnover, T s — T h , versus B, for all magnetars for which both T s 
and Th are measured. Circles represent pulsed flux, triangles represent total flux. See text 
and Table 1 for data used and references. Right: Same but versus v. Both: Error bars 
represent la uncertainties. All known RPPs for which the quantity has been measured have 



T s - T h « 0, with B in the range 1 - 10 x 10 12 G. 



